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We demonstrate that the time-dependent Krieger-Li-Iafrate approximation in combination with 
the exchange-only functional violates the 'Zero-Force Theorem'. By analyzing the time-dependent 
dipole moment of Nas and Na^j", we furthermore show that this can lead to an unphysical self- 
excitation of the system depending on the system properties and the excitation strength. Ana- 
lytical aspects, especially the connection between the 'Zero-Force Theorem' and the 'Generalized- 
Translation Invariance' of the potential, are discussed. 
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Since the pioneering work of Ando, Peuckert, Zang- 
will and Soven [J, [3, S| find the rigorous foundation by 
Runge and Gross time-dependent density functional 
theory (TDDFT) has become one of the most success- 
ful theories for the calculation of electronic properties. 
Based on the idea of describing a system only in terms 
of its time-dependent density, TDDFT offers the oppor- 
tunity to access large systems which are out of reach for 
any wave-function based method. Additionally, it offers a 
tractable method to deal with strong, non-linear and non- 
perturbative excitations of atoms, molecules and clus- 
ters. High-harmonic generation and above-threshold ion- 
ization are just two examples for applications of TDDFT 
in this domain. 

As in ground-state density functional theory, most 
TDDFT calculations are done in the Kohn-Sham (KS) 
scheme In this approach, the interacting many- 

particle system is replaced by a ficticious system of non- 
interacting particles moving in a local effective potential 
wks- This potential is chosen in such a way that it ex- 
actly reproduces the time-dependent density. As in the 
static case, the potential is split into the Hartree part 
wh, containing the classical electrostatic interaction, the 
local external potential Woxt from the ions and other ex- 
ternal potentials, e.g., from a laser, and the exchange- 
correlation (xc) potential fxc containing all quantum- 
mechanical exchange-correlation effects. Since the ex- 
act xc potential is not known, it is crucial to have good 
approximations for it. The most well-known approxima- 
tion is the time-dependent local-density approximation 
(TDLDA). It simply uses the static local-density approx- 
imation functional for v^c in combination with the time- 
dependent density. Although being based on the ground- 
state energy of the homogenous electron gas, this approx- 
imation works well in many situations [s, [a, S B S ■ 



Nevertheless it can also fail dramatically in some situa- 
tions (see, e.g., U, 12, One of its most prominent 



problems is the self-interaction error leading to a wrong 
asymptotic behavior of Uxc and, as a consequence, wrong 
ionization dynamics. In addition, it lacks a derivative 
discontinuity and does not include any memory 

effects 1^. 

One possible way of overcoming these problems are 
orbital functionals 17]. These functionals depend ex- 
plicitly on the KS orbitals and implicitly on the den- 
sity. Thus, they are legitimate density functionals and 
provide many advantages over explicit functionals of the 
density. For instance, they include memory effects, show 
a derivative discontinuity, and using the exact-exchange 
functional (EXX) with the KS orbitals cures the Hartree 
self-interaction problem. Unfortunately, they also pose 
a severe problem: in order to find Wxc, it is necessary 
to solve the time-dependent optimized-effective poten- 
tial (TDOEP) equation [l3|. This is not an easy task 
[igt - To circumvent this problem, Ullrich et al, already 
in the first publication of the TDOEP equation, pro- 
posed an approximation to the exact potential, namely 
what they called the time-dependent Krieger-Li-Iafrate 
(TDKLI) approximation. TDKLI has been frequently 
used to carry out TDDFT calculations with orbital func- 
tionals, e.g., to calculate ionization processes and high- 
harmonic generation 

[11 S EH- In fact, besides the 
time-dependent common energy denominator approxi- 
mation (CEDA) 13, 2^, this is the only practicable ap- 
proximation to the exact TDOEP outside the linear do- 
main at present. On first sight it might be surprising 
that although the EXX-TDKLI potential is not obtained 
as the functional derivative of the EXX orbital functional, 
it is nevertheless a legitimate density functional due to 
the one-to-one correspondence between the density and 
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the external potential applied to the non-interacting KS 
system. 

On the fundamental side of TDDFT, many exact 
constraints which v^c must fulfill have been revealed 
24, [2^ [2^. Well-known examples are the 'Harmonic- 
Potential theorem' 27| and the 'Zero-Force theorem' 
@i 28, 2§|. The statement of the latter is that the force 
which the xc potential exerts on an electron cloud must 
vanish, i.e., 



n(r,t) Vt;xc(r,0 d^r 







(1) 



Since the force from the Hartree potential vanishes, this 
guarantees that 



n(r,i) Vwoxt(r,i) d^r 



holds with P{t) being the electronic momentum of the 
system |25| ■ From earlier work it is known that the EXX 
functional is a 'conserving approximation' [26| if the full 
TDOEP equation is solved. However, since this is a cru- 
cial requirement for the proof, the situation is less clear 
for the case of the EXX-TDKLI potential. It is the aim 
of this manuscript to explore whether the EXX-TDKLI 
nevertheless satisfies the 'Zero-Force theorem'. 

In order to answer this question, we have calculated the 
response of a Nas cluster to a small dipole excitation. 
This cluster has planar geometry (oriented here in the 
{x — j/)-plane) and a rather "soft" electron cloud which 
provides a critical test case for our purposes. For the 
calculation we used a modified version of the PARSEC 
program [30| in combination with a local pseudopotential 
for the Na cores 3l|. For the real-space grid a sphere of 
radius 20 uq and a grid spacing of 0.7 oq was used. The 
kinetic energy part of the KS Hamiltonian was approxi- 
mated by a 12-point formula for the Laplacian. After the 
ground-state calculation, we applied a momentum boost 



exp(ir -pboost/Ti) with |pbo 



3.834 X 10"*;ian 



responding to a total excitation energy of the system of 
1.0 X 10~^ eV, to all KS orbitals. The boost had an equal 
strength in x-, y-, and ^-direction. The resulting excited 
state was propagated in real time with fixed ions. The 
propagation was done with a fourth-order Taylor approx- 
imation to the propagator. We checked the convergence 
of our results with respect to the time step and found a 
value of At = 0.003 fs for the propagation sufficient. 

Fig. [1] shows the resulting y-component of the dipole 
moment dy(t) = e Jynir^t) (Pr (with e being the elec- 
tron's charge). For a while the amplitude shows reason- 
able oscillations in agreement with the initial boost. But 
after about 400 a.u., it increases rapidly and steadily. For 
comparison, the same quantity is plotted for a TDLDA 
calculation. There, no increasing amplitude is observed. 
Since the TDLDA satisfies the 'Zero- Force theorem' [29| . 
this is already a hint that the increasing amplitude in 
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FIG. 1: (Color online) y-component of the dipole moment of 



(2) Nas after an initial boost of |pboost| = 3.834 X 10"* hag 



In contrast to the TDLDA result, the EXX-TDKLI curve 
shows an increasing amplitude of the oscillation. The EXX- 
TDKLI curve is slightly shifted for better comparison. Ryd- 
berg atomic units are used throughout. 
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FIG. 2: (Color online) Time derivative of the y-component of 
the total momentum of the electron density. 'Woxt + «x' and 
'iicxt' label the results calculated via Eq. ^ with the exter- 
nal potential or, respectively, the external and the exchange 
potential. The violation of the 'Zero-Force theorem' is visible 
in the difference between 'iiext + fx' and 'noxt'- 



the TDKLI calculation is related to a violation of the 
'Zero- Force theorem'. 

To demonstrate that this is indeed the case, we have 
monitored the expected time-derivative of the total mo- 
mentum, i.e., we have calculated the right-hand side of 
Eq. ^ with the external potential coming from the local 
pseudopotentials of the ions [31]. 

In addition, we have evaluated the left-hand side of 
Eq. ([T|), which should be zero. Fig. [5] shows the result 
obtained from Eq. ^ and the sum of Eq. and Eq. 
([T]). One observes that the total force from the external 
potential on the electron density, Eq. differs signifi- 
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cantly from the total force obtained from the sum of both 
equations. This clearly demonstrates that Ux violates the 
'Zero-Force theorem' and contributes to the total force. 
Since the time-dependent dipole moment is connected to 
the total force F{t) via (see, e.g., (25j l 



d^d{t) = e9tP(<) = eF(t) 



(3) 



it follows that the violation of the 'Zero-Force theorem' 
leads to a wrong time-dependent dipole moment and to 
a self-excitation of the system. Finally, we have checked 
that Eq. ^ holds for the KS system. For this we calcu- 
lated the force from the total KS potential, i.e.. 



F{t) := - J n(r,i) VwKs(r,i) 



(4) 



and compared it to the time-derivative of the KS current 



jKs(r, t) 



^2 N 
k=l 



c.c.) (5) 



from our calculation. The current is connected to the to- 
tal momentum via P(t) — m J jKs(r, t) (Pr. By explicitly 
calculating the second time-derivative of the monitored 
dipole signal, we confirmed that Eq. ([3]) holds in our cal- 
culation as it should be. 

To study the influence of the 'Zero- Force theorem' vi- 
olation in more detail, we have carried out the same 
calculation for a Naj]" cluster. In contrast to Nas, this 
cluster can be considered as one of the most 'forgiving' 
systems because of its spherical shape and the positive 
charge leading to a stable 'plasmon'-like oscillation when 
excited. Fig. [3]shows the resulting violation of Eq. A 
close look reveals that the 'Zero-Force theorem' is slightly 
violated again. But, in contrast to Nas, the violation 
does not increase in time. Checking the time-dependent 
dipole moment shows that it is also stable. Beside the 
stronger binding forces a possible explanation for this 
observation could be that the higher inversion symmetry 
of the almost spherical Na^ leads to an error cancella- 
tion in the course of one density oscillation and thus to 
a strongly reduced increase of the violation. In any case, 
the result clearly corroborates the intuitive expectation 
that the system properties have a strong influence on the 
degree of the violation of the 'Zero-Force theorem'. 

In addition to the system properties, one can also ex- 
pect an influence of the excitation energy on the 'Zero- 
Force' violation. And indeed, this can be found. Fig. 
m shows the results for Na5 excited with three different 
boost strengths of 3.834 x IQ-'^Tia'^^, 6.062 x IQ-'^fia^^, 
and 8.573 x 10^'*?iaQ^. Obviously, the deviation from 
zero varies with the boost strength. The importance of 
the excitation energy can also be demonstrated by con- 
sidering the extreme situation of no excitation at all. In 
this case, the time-dependent KS orbitals are given by 
the ground-state orbitals multiplied by a time-dependent 




100 200 300 400 500 600 
Time [a.u.] 

FIG. 3: (Color online) Violation of the 'Zero-Force theorem' 
for Nag" and Nas. Again, the j/-component is plotted. Both 
systems were excited by a momentum boost corresponding to 
an excitation of the system by 1.0 x 10~^ eV. In contrast to 
Nas, the resulting curve for Nag" does not increase in time. 
The resulting dipole moment of Na^ is also observed to be 
stable. 
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FIG. 4: (Color online) Violation of the 'Zero-Force theorem' 
for 

6.062 X 10" 



Nas excited with different boosts of 3.834 x 10 *' ha^^ , 



ha^ , and 8.573 x 10~ ha^ . Higher excitation 
energies clearly lead to an increased violation. As in Fig. [2] 
only the j/-component is plotted. 



phase factor exjp(—iekt/h) containing the KS eigenvalue 
ek- But since this phase factor does not influence the 
potential, the whole system remains in a stationary state 
and the violation of the 'Zero- Force theorem' remains 
constant. This constant violation does not lead to a 
non-stationary state since / n(r, <:) VwKs(r, t) (fir van- 
ishes due to the ground-state iteration. In other words, 
for the ground-state density the total force from the ex- 
ternal and the xc potential are in equilibrium leading to 
a stationary state. 

Beside the just discussed aspects, there are three other 
conditions under which a violation of the 'Zero-Force 
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theorem' in a TDKLI calculation may not be observed. 
First, for short time scales the accumulation of the vi- 
olation can be too small to show up significantly. Cer- 
tainly, this time scale depends on the two aspects dis- 
cussed above: the excitation strength and the system 
properties. The second situation occurs when a strong, 
ionizing external field is applied to the system. This can 
hide the error in the xc potential completely. Finally, 
for spin-saturated two-particle systems the TDKLI and 
the TDOEP potential coincide and, as a consequence, 
the 'Zero- Force theorem' is rigorously satisfied if the TD- 
KLI potential comes from a 'conserving approximation' 



32] ■ Since, to the best of our knowledge, one of these 
three conditions can be found in all applications of the 
TDKLI potential to date, it is not surprising that the 
violation of the 'Zero-Force theorem' was not reported 
earlier. 

After the numerical results we now focus on analyti- 
cal considerations. As Vignale has shown, any potential 
obtained from a 'generalized-translation invariant' xc ac- 
tion functional satisfies the 'Zero- Force theorem' 2^, 2^. 
In addition, the potential also satisfies 'Generalized- 
Translation invariance', i.e., v^c rigidly follows a rigidly 
translated density, and, as a consequence, the 'Harmonic- 
Potential theorem'. Actually, the same arguments show 
that if Vxc satisfies 'Generalized- Translation invariance' 
and is the functional derivative of some xc action func- 
tional, it must also satisfy the 'Zero-Force theorem'. 
Since the EXX-TDKLI satisfies the 'Harmonic-Potential 
theorem' [^, [s^, we conclude from our numerical re- 
sults that the EXX-TDKLI potential in general cannot 
be obtained as the functional derivative of some xc ac- 
tion functional with respect to the density. This result is 
in line with earlier results for the static Slater potential 
SJI . The missing action functional for the EXX-TDKLI 
potential, in combination with the observation that the 
TDKLI potential is nevertheless used as part of a legit- 
imate density functional procedure, makes any rigorous 
analytical statements about the properties of the poten- 
tial extremely difficult. A hand- waving argument, how- 
ever, is provided by the expression for the EXX-TDOEP 
involving the orbitals and orbital shifts [3, Since 
the full expression satisfies the 'Zero- Force theorem', it 
is not surprising that the TDKLI approximation, which 
neglects the terms containing the orbital shifts, violates 
the 'Zero-Force theorem'. Additionally, the violation of 
the 'Zero-Force theorem' in the static Slater approxima- 
tion [s^ supports the presented numerical results, too. 

Finally, we want to comment on the generalization of 
our results to TDKLI potentials from other xc orbital 
functionals for the action. Clearly, the possibility that 
other expressions lead to a TDKLI potential which sat- 
isfies the 'Zero- Force theorem' cannot be ruled out from 
our calculations. However, since the EXX functional is 
a 'conserving approximation' [26| and this property is 
destroyed by using the TDKLI and not the TDOEP po- 



tential, it is highly plausible that this also happens if 
other 'conserving approximations' (26l . [s^ l are used. In 
any case, it is highly recommendable to check the 'Zero- 
Force theorem' in any TDKLI calculation. 

To summarize, we have demonstrated that the TDKLI 
approximation in combination with the exact-exchange 
functional violates the 'Zero- Force theorem'. By calcu- 
lating the response of Nas and Nag , we have demon- 
strated that this violation can lead to a self-excitation of 
the system and a wrong time-dependent dipole moment. 
Furthermore, a strong system and excitation-energy de- 
pendence has been observed. Due to the difficulties as- 
sociated with the construction of orbital-dependent ap- 
proximations for the xc potential which satisfy the 'Zero- 
Force theorem', our findings clearly demonstrate the ur- 
gent need for a working TDOEP scheme in order to make 
orbital functionals accessible for real-time TDDFT. 
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